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A phase separation model is presented for the stripe phase 
of the cuprates, which allows the doping dependence of the 
photoemission spectra to be calculated. The idealized limit 
of a well-ordered array of magnetic and charged stripes is 
analyzed, including effects of long-range Coulomb repulsion. 
Remarkably, down to the limit of two-cell wide stripes, the 
dispersion can be interpreted as essentially a superposition 
of the two end-phase dispersions, with superposed minigaps 
associated with the lattice periodicity. The largest minigap 
falls near the Fermi level; it can be enhanced by proximity to a 
(bulk) Van Hove singularity. The calculated spectra are dom- 
inated by two features - this charge stripe minigap plus the 
magnetic stripe Hubbard gap. There is a strong correlation 
between these two features and the experimental photoemis- 
sion results of a two-peak dispersion in La2-a:Sra;Cu04, and 
the peak-dip- hump spectra in Bi2Sr2CaCu2 0g+i . The differ- 
ences are suggestive of the role of increasing stripe fluctua- 
tions. The 1/8 anomaly is associated with a quantum critical 
point, here expressed as a percolation- like crossover. A model 
is proposed for the limiting minority magnetic phase as an iso- 
lated two-leg ladder. 



I. INTRODUCTION 

Evidence for stripe phases in the cuprates continues 
to grow. Particularly in the La2_xSr2;Cu04 (LSCO) 
family, a convincing case for (predominantly dynamic 
or disordered) stripes can be flpftde, based on ela^ 
and inelastic neutron scatteringSiS, NMR and NQI 
In other systems, the evidence is moK|-|ambiguous. In 
YBa2Cu307_5 (YBCO), there is nowBa clear evidence 
for incommensurate modulation of the inelastic mag- 
netic neutron scattering near Q — (tTjTt), but so far onijt 
in underdoped YBa2Cu306.6- Balatsky and Bourgesty 
find a broad commensurate peak, but the width of the 
peak scales with doping in exactly the same way as 
the incommensurability in LSCO, suggestive of an un- 
resolved underlying incopimensurability in YBCO, as 
well. Also, de LozanneO finds direct STM evidence 
for incommensurate modulations (parallel to the chains), 
with a similar periodicity to the neutron data. MookO 
has reported similar incommensurate neutron peaks in 
BiaSrzCaCuaOg+i (BSCCO). Potentially stripe-related, 
phonon anomalies have been reported in both LSC0t2l 
and YBCdMn Doping with Zn seems ±a|-S|tabilize the 
stripe phaseE-3. Photoemission evidenceE3ll3 for stripes 
has been controversialMltj. 

Over the same doping regime, there is also evidence 
for a pseudogap, and it is an important problem to un- 
derstand how both pseudogap and stripes can coexist. 



In particular, photoemission finds a dispersion consistent 
with the two-dimensional (2d) energy bands, whereas in 
the stripe phase the magnetic stripes should be insulat- 
ing, leading to a one-dimensional (Id) dispersion along 
the charged stripes. 

The presence of stripe phases raises important issues 
of how energy dispersion and even Fermi surfaces can 
be well defined concepts in the presence of fluctuating 
stripes. An important insight li^to this problem is the 
finding by Salkola, et al. (SEK)t3 that a well-defined av- 
erage dispersion persists even in the presence of strongly 
fluctuating stripe orider. However, in that paper, and 
a related calculationEll, the stripes were modelled by a 
charge density wave like order, with sinusoidallyjiaQing 
hole density Several unrestricted Hartree-FockEScJ or 
slave bosontj calculations find evidence for much sharper 
density variations. The present paper analyzes a phase 
separation scenario, modelling the stripes as associated 
with free energy minima at two characteristic hole den- 
sities. This allows the doping dependence of the stripes 
and the resulting photoemission spectrum to be analyzed. 

It is found that long-range stripe order can persist 
even in the presence of Coulomb interactions. The re- 
sulting dispersion is clearly recognizable as a superposi- 
tion of the magnetic and charged stripe dispersions, with 
superimposed minigaps due to the stripe order. These 
dual dispersions provide a natural interpretation for the 
experimentally observed photoemission dispersions, ty- 
ing together results on LSCO, BSCCO, and SraCuOaCla 
(SCOC). In the model, the 1/8 anomaly can be under- 
stood as a form of quantum critical point (QCP), asso- 
ciated with a crossover between a magnetic stripe dom- 
inated regime and a charged stripe dominated regime. 
Within the latter regime, the (tt, tt) spin gap in YBCO 
is related to the behavior of a two-leg ladder (isolated 
magnetic stripe). 

Remarkably, within the charge stripe dispersion, a 
clear signature of the two-dimensional Van Hove singular- 
ity (VHS) persists, down to the limit of a single, two-Cu 
wide stripe. There is a strong coupling of the minigaps 
with this VHS, leading to a novel stripe-induced VHS 
splitting. The doping dependence of this splitting closely 
resembles that of the pseudogap. 

The paper is organized as follows. Section II shows 
that a low hole doping of the charged stripes, xq ^ 0.25 
is not only compatible with experiment, but also makes 
sense theoretically, in terms of kinetic-energy stabilized 
stripes. The models for the magnetic and the charged 
stripes are introduced in Section III, along with a dis- 
cussion of long-range Coulomb interaction. Section IV 
gives the results of the stripe calculations, which self- 
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consistently determine the hole distribution. The dop- 
ing dependence of the dispersion is presented, for vary- 
ing strengths of Coulomb repulsion. Finally, the effect 
of an additional (ferromagnetic) interaction on splitting 
the VHS degeneracy on the charged stripes is discussed. 
In Section V, these results are compared to experiment, 
and a consistent model of the photoemission in LSCO 
and BSCCO is presented. Section VI points out that the 
model has a QCP - actually, a series of 'magic doping' 
QCP's, of which the simplest is the 1/8 anomaly. To il- 
lustrate the resulting crossover, an additional calculation 
is presented in Section VII, showing how the doping de- 
pendence of the spin gap in YBCO can be understood. 
Possible explanations are also presented for the satura- 
tion of the incommensurability S vs x found by Yamada, 
et al.H 

Section VIII includes discussions of the interpreta- 
tion of the pcak-dip-hump structure in BSCCO, the new 
stripe- VHS spitting pseudogap, a discussion of Fermi sur- 
faces and remnant Fermi surfaces in the stripe phase, and 
comparison with earlier calculations. A summary of the 
principal conclusions of this work is given in Section IX. 

II. FRACTIONALLY-OCCUPIED STRIPES 

A. Comparisons with other Oxides 

Stripe arrays have now been found in a number of 
oxides, most notably nickelates and manganites. The 
similarities of cuprates with nickelates are particularly 
close: in both systems, the charged stripes act as an- 
tiphase boundaries for the magnetic stripes, and in both, 
the charge ordetj-arises at higher temperature than the 
magnetic ordernnZl. The nickelate stripes run diago- 
nally (with respect to the Ni-O-Ni bonds); this is also 
true of the LSCO stripes, in the spin glass regimec3, 
X ~ 0.04—0.06. However, in the superconducting regime, 
X > 0.06, the cuprate stripes are generally horizontal and 
vertical. 

One striking difference is that in the nickelates and 
manganites, the charged stripes correspond to integer 
doping . (one hole per Mn or Ni), leading to simple 
patternsE^cj of commensurate stripe arrays. There are 
prominent phase transitions at rational fractions, 1/2, 
and 1/3, corresponding to holes on every nth row, with 
evidence for commensurability locking in between (i.e., 
the 1 /3 phase persists in an extended doping range about 
X = 1/3.) Consistent with integer filling, the phases are 
all insulatin^il. In contrast, in the cuprates the phases 
are all conducting or weakly localized, and the only frac- 
tion which appears prominently is 1/8. 

In the present paper, a simple explanation is proposed 
for this distinction. The charged stripes are fractionally 
doped, with approximately 1/4 hole per Cu (hence ex- 
plaining the finite conductivity) . The magical 1 /8 doping 



would then correspond to the simplest 'commensurate' 
pattern of these stripes. 

The stability of the stripe phase decreases in the se- 
quence manganites, nickelates, cuprates. Thus, while 
there are beautiful electron microscopic, images of long- 
range stripe order in the manganiteso, stripes in the 
cuprates are mainly fluctuating, with only short-range 
order. Within the present model, this pattern is read- 
ily understood, since the charged stripes are stabilized 
by CDW instabilities; this-is similar to models for the 
nickelates and manganitest3. The strength of this in- 
stability can be estimated by comparing the strength 
of electron phonon coupling, which follows the same se- 
quence: manganites (with well-defined Jahn- Teller po- 
larons), nickelatesu, cuprates. It is only in the cuprates 
where the interaction is so weak that a fractional occupa- 
tion can be stabilized, and it is only in the cuprates that 
the stripe formation is so weak that superconductivity 
can successfully compete. 

B. Origin of Fractional Occupation 

Hartree-Fock calculationsEl of the t J model find that 
the holes condense onto domain walls between antifer- 
romagnetically ordered domains, producing fully occu- 
pied charge stripes - one hole per Cu. However, neu- 
tron diffractionEJ finds a charge modulation of periodicity 
four Cu atoms a,t x = 0.125_which implies only 1/2 hole 
per cell. Tranquada, et al.El suggested a model for the 
charged stripes, based on their experience with stripes 
in nickelates. The hole-doped stripes are one cell wide, 
and have a hole on every other site. A rpicroscopic model 
for such a domain wall can be derivedEa by incorporat- 
ing a charge-density wave (CDW) instability along the 
stripes, treating them as one-dimensional metals. How- 
ever, such states with integral hole doping are likely to 
be insulators, as is the case in the stripe phases of the 
nickelatestil, whereas the cuprates are either conducting 
or weakly localized. 

Moreover, fractional hole occupation would seem to be 
more natural for the tJ and Hubbard models, since the 
energy of doped holes is lowered, by finite hopaing t in 
a partially filled band. Visschert^ and NagaevE3 showed 
that the holes enhance their kinetic energy by creating 
local ferromagnetic domains (ferrous) in which they are 
free to hop. This leads to a preferred hole density, Xf in- 
side the ferron domain. In a two-dimensional, tJ version 
of the model (letting Ti^ /2m to? ^ with a the lattice 
constant), 

= 0.334, (1) 

with z — A the number of nearest neighbors of a given 
Cu, and I have assumed J/t ^ P-?5. A similar result 
was found by Nayak and Wilczektll. Nagaav's model is 
a large-S* theory, and Emery and Kivelsoi£a extended it 
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to 5* = 1/2, although thes2-.did not address the issue of 
X f . Auerbach and LarsonEj showed that a single doped 
S — 1/2 hole will spread out over a ferromagnetic domain 
covering 5 lattice sites, suggesting a comparable value 
for X/, > 0.2 holes per cite, on average (since the hole 
has a higher probability of being on the central atom). 
Recent density matrix renormalization group (DMRG) 
calculations of th£ tj model undertaken by White and 
Scalapino (WS)e2IcJ find charged stripes which are two 
Cu's wide, with an average hole doping of 0.25 hole per 
Cu on the charged stripes. These calculations are further 
discussed in Appendix A. 

Since the charged domains are stabilized by the hole 
kinetic energy, it is plausible that enhancing the ki- 
netic energy could enhance the stability of the hole- 
doped stripes. Thus, in a generalized Hubbard model, 
with next-nearest neighbor hopping t' , it is found that 
a macroscopic ferromagnetic phase is stabilized in the 
vicinity of the Van Hove Singularity|-(iVHS)c3. Moreover, 
an extended Hartree-Fock analysiscJ finds phase sepa- 
rated states smoothly evolving between the AFM and FM 
regions, from a single magnetic polaron to FM stripes to 
a uniform FM phase. 

However, such ferromagnetic domains have not been 
observed in the cuprates. Nevertheless, there are al- 
ternative VHS routes to fractionally-occupied stripes. 
The large density of states (dos) associated with a 
VHS can d rive . a large number of competing electronic 
instabilitiescScJ, and it was early suggested that this 
could be thje-| origin of nanoscale phase separation in 
the cuprateala. In particular, it was demonstrated that 
strong electron-phonon coupling cpuM stabilize a charge- 
density wave phase near the VHScJ'cS. 



physics of the 2-leg ladder. When one eliminatescS cer- 
tain one-dimensional operators (which break the fc — > — fc 
symmetry along the ladder), one is left with the same 
S0(6) group introduced earlier for the VHS. Such a cor- 
respondence would fail for a single-leg ladder. 

Hence, the present model is restricted to stripes 
which are an even number of cells (or Cu atoms) wide. 
This point was previously postulated for the magnetic 
stripes,|-iii terms of spin gaps associated with even-legged 
laddersEj. Moreover, WS find tss40-Cu wide charged 
stripes in their DMRG calculationsEj. With this assump- 
tion, it is found that a VHS-like feature can be clearly 
resolved near the Fermi level in the stripe phases. More- 
over, the stripes provide a new mechanism for VHS split- 
ting - minigaps - which can generate a pseudogap with 
the correct doping dependence. 

III. MODELING THE STRIPES 

While the stripes are likely to be strongly fluctuating, 
the band structure modifications should be strongest, 
and can be analyzed in most detail, in an ordered stripe 
phase. Hence, the present calculation assumes perfectly 
ordered stripe phases to describe this 'worst case' sce- 
nario. It will be assumed that there are two preferred 
hole densities, a; ~ on the magnetic stripes, and xq ~ 
0.25 holes per Cu on the hole-doped stripes. Coulomb 
effects lead to additional charge relaxation, and a more 
uniform distribution of charge. Section III.D. 

A. Model for the Magnetic Stripes 



C. Viability of VHS Models 

In any model of stripe phase formation based on Fermi 
surface features, there is a fundamental question of self- 
consistency: do the features persist in the limit of an iso- 
lated stripe? Can one still recognize bulk features of the 
band structure and Fermi surfaces of the phases forming 
the stripe array? This is one of the main issues that this 
paper resolves: even in the limit of nanoscopic stripes, 
the band structure is recognizably a superposition of the 
structures of the two end phases. The main role of stripe 
order is to introduce miniband gaps into this structure. 

In the particular case of the VHS's, there were a num- 
ber of preliminary indications jSSihich suggested such an 
affirmative answer. First, SEKE3 found that an average 
dispersion persists in the presence of fluctuating stripes; 
the resulting 'flat bands' are a signature of the VHS, Sec- 
ondly, within a group theoretical (SO (6)) modelcj, the 
Van Hove instabilities all remain well-defined on a sin- 
gle plaquette of 2x2 Cu atoms, so a fortiori they should 
remain well defiuad on a 2-leg ladder. Indeed, Lin, Ba- 
lents, and FishciO found an SO (8) group controlling the 



In the insulating phase, a variant of the spin- 
density wfiY\ (SDW) model studied by Schricffer and 



coworkersElH is used. |— .This model works surprisingly 



well in the large-U limitE2l, reproduces the spin wave spec- 
trum of the Heisenberg model, and has served as the 
basis for a_number of extended treatments of correlation 
effectsCl^ES. For realistic parameters {t, t' , U), the model 



has a Mott-Hubbard gap of 2eV, a,a 
dispersion found in the oxycloride; 



i^an reproduce the 
Appendix B. 



The dispersion of the one-band model can be written 



At CxCy 



(2) 



with Ci = cosfcifl. Writing e± — (efc±efc+Q)/2, the eigen- 
values in the presence of a Hubbard U become 



± J el + C/2, 



(3) 



where U ~ Unig. In the limit U >> t, the lower Hub- 
bard band may be approximated 



-U- 



At Cj-Cy 



J{Cx 



(4) 
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with J = 2t^ /U . The parameters can be deter- 
mined by fitting to the observed photoemission disper- 
sion in SCOC. For simphcity, one can use analytical 
expressions for the parameters at three fc-space points: 
^-(7r/2 7r/2) ^ -^ U , E-{tt,0) ^ -~U + At' , E~{0,0) = 
—At' — \/U^ + I6t^ (Eq. 1^ is not sufficiently accurate for 
this purpose). The fit yields t — 325meV, U/t — 2.5, 
and T = 2t'/t — —0.552. Solving the gap equation at 
half filling, this value of U corresponds to U/t — 6.03, 
Mq{x = 0) = 0.414 (Fig. |l|), or 83% of the classical 
value. 

For these parameters, M{x) is multivalued for x > 
0.38. This implies that the magnetic to non-magnetic 
transition is first order. This is discussed further in Ap- 
pendix B. However, this density is rather higher than ex- 
pected for charged stripes. In LSCQ, the VHS splitting 
seems to terminate near x = 0.26E3E3, and similar resailts 
are found below for YBCO. In a number of modelsCJO, 
the AFM instability is replaced by a second instability, 
driven by splitting the VHS degeneracy. Note that the 
bare {U = 0) VHS falls at x = 0.25 for r = -0.559, very 
close to the value needed to explain the dispersion in the 
insulating phase. 



M 




B. Model for the Charged Stripes 

It is assumed that the hole-doped stripes are stabilized 
by splitting the VHS degeneracy, at the doping xq ~ 0.25 
where the VHS falls at the Fermi level. An earlier slave 
boson calculationc3 demonstrated that electron-phonon 
coupling could provide that stabilization energy, even in 
the presence of |Strong correlation effects. A ferromag- 
netic interactioro can produce similar splitting. 

While the earlier electron-phonon calculation involved 
a three-band model, here a simpler one-band model will 
be adopted. A parametrized form of the freei-energy vs 
doping found in the self-consistent calculation^, Fig. ||, 
will be assumed, to stabilize the stripe phase. It is conve- 
nient at present to not introduce any mechanism to split 
the VHS degeneracy. This allows a definitive answer to 
an important question: can evidence for the VHS still 
be found in the presence of a well-defined stripe phase? 
The answer is a clear yes: the resulting dispersion is a 
superposition of the magnetic dispersion and the charged 
stripe dispersion, with recognizable VHS feature. What 
is more, the stripe phase minigaps provide a new mech- 
anism of VHS splitting, with a doping dependence com- 
parable to the experimental pseudogap. 

A very simple doping dependence of the parameters is 



assumed. From Eqs. 
a factor t/Umq, Fig. 



i, for finite U t is renormalized by 
iT so the increase of t with doping 
is accomplished by the decrease in niQ, the ordered mo- 
ment. We will thus make a simple Ansatz that the only 
effect of doping is to renormalize 



mg toq(1 - x/xq). 



(5) 



FIG. 1. Doping dependence of magnetization in the SDW 
model. 



Since the stripes are predominantly near the limiting 
states X = 0, xq, the detailed nature of the intermediate 
states is relatively unimportant. As noted above, Eq. ^ 
neglects the gap on the charged stripe; in Section IV. D, a 
ferromagnetic interaction will be included on the charged 
stripes, to show that the VHS splitting is preserved in the 
striped phase. 



If Mq is interpreted as the long-range antiferromag- 
netic order parameter, then the model does a poor job 
in describing the temperature and doping dependence of 
the Neel transition^, Tn, yielding Tn ~ C//4. Figure |l| 
shows that, while Mq is strongly renormalized by doping, 
the mean-field theory underestimates the rapidity of the 
falloff of Tpf with x. However, the mean-field results are 
best reinterpreted as representing short-range order - the 
magnetic fluctuations ~ and hence the renormalization of 
the splitting into upper and lower Hubbard bands. In 
this case, the mean- field calculations are injgpod agree- 
ment with exact diagonalization calculationsEj. The fact 
that the gap is much smaller in the dopei phase is consis- 
tent with the experimental observationtJ that the upper 
Hubbard band rapidly disappears with doping. 



C. Free Energy Minima 

To stabilize the stripe densities at the values a; = 
for magnetic stripes, and x = xq = 0.25 for the charged 
stripes, the following free energy is introduced, based on 
the results of slave boson calculations for the three-band 
modelcj: 

fo{x) = ^lox{^ - -)^ (6) 

Xo 

for X > 0. (At a; there is a cusp-like minimum, 
associated with the Mott gap in the chemical poten- 
tial.) Figure ^ illustrates the free energies calculated 
in a three-band slave boson calculation for competing 
magnetic (here a flux phase - dashed line) and charged 
(CDW - solid line) phases, and a fit of these to Eq. ||, 
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with fiQ = 0.9eV. For these calculations, a;o = 0.16 was 
assumed. This should be compared to Fig. 2 of Ref. 
Note that there is an error in the caption of that figure: 
the CDW phase there corresponds to a weaker coupling, 
Vep = 0.6eV. 



f(eV) 




FIG. 2. Jjree energy calculated in a three-band slave bo- 
son theoryCj, showing phase separation between a flux phase 
(dashed line) and a CDW phase, with Vep ~ leV (solid line). 
Dotted line is a fit to Eq. n 

Equation ^ is a convenient form for parametrizing the 
confining potential of a striped phase. It has only two 
parameters, xq and /iq, or equivalently, = 4/Ltoa^o/27, 
the maximum free energy barrier, at xq/S. In the present 
stripe phase calculations, these parameters are taken as 
xo = 0.25 and = 0.312ey, or /,„ — llmeV. This value 
corresponds to Vep — 0.6eV of Ref. |4^, and allows us to 
see that even a relatively modest confining potential can 
stabilize the stripe phase against the Coulomb potential. 
This free energy corresponds to an additional chemical 
potential 



Xo 



-)(1 



3a; , 

Xq' 



(7) 



for X > 0. In the calculations, this fi{x) is added to 
the potential on each row, and the local density adjusted 
until self-consistency is attained. 

At X ~ 0, fi has a discontinuity, the Mott-Hubbard 
gap. Hence, at this point, the Fermi level can take on any 
value inside the gap. To model this in a computationally 
stable manner, the discontinuous step in fi is replaced by 
a linear ramp, connecting the values of /i at x = —0.01 
and X = -1-0.01, and assuming ^{x = —0.01) = — 
-hO.Ol). Thus, when the calculation finds \x\ < 0.01, 
it generally implies that the Fermi level is in the gap 
of the magnetic stripes. However, due to hybridization 
with holes in the charged stripes, it is possible to have a 
well-defined Mott gap, with a small doping x > on the 
magnetic stripes (typically, x < 0.05). 



D. Madelung Energies of Stripes 

We will assume for simplicity that all stripes, both 
magnetic and charged, are an even number of cells wide. 
This means that only a relatively small number of stripe 
configurations are involved in the doping range of inter- 
est. For instance, labelling the stripe configuration by 
m, n, where m is the width of a magnetic stripe and n 
the width of a charged stripe, we will explore in detail the 
pure phases m,n = 6,2 {x = xo/i ~ 0.0625, if xo ~ 0.25), 
4,2 (x = xo/3 ~ 0.0833), 2,2 (x = xo/2 ~ 0.125 
- the 1/8 phase), 2,4 (x = 2xo/3 ~ 0.167), and 2,6 
(x = 3xo/4 ~ 0.1875). Intermediate dopings would cor- 
respond to mixed phases. For each of these phases, we 
assume that there can be different dopings on each row; 
by symmetry, there can be (m + n)/2 inequivalent rows 
for the m,n-phase. 

In the presence of charging, it is the electrochemical 
potential fie and not the chemical potential, /i, which is 
constant. For electrons, fie — fJ- — eV, where V is the 
electrical potential. Given the average hole density on 
each row, V can be calculated as a Madelung sum. For 
each configuration, the Madelung sum can be calculated 
for each row. Actually, since the overall chemical poten- 
tial must be adjusted to fix the total hole density, all that 
need be calculated is the difference in Madelung potential 
between the different rows. This is calculated by assum- 
ing a pure Coulomb interaction, screened by a static di- 
electric constant, e. The on-site term is neglected, having 
already been included as U. 

The various Madelung constants can be expressed as 
follows. For the (6,2) stripe, label the rows 1,2,3,4, with 
4 = the charged stripe, and 1 (3) = the magnetic rows 
farthest from (nearest to) row 4. Let Vi be the Madelung 
potential for the ith row, Vi = (Vi — Vi)xo, Xi the hole 
doping of the zth row, and Xi — (x^ — Xi)/xo. Then 



(8) 



where the X'" matrices have been calculated numerically, 
with results listed in Table I, for the cases (m,2), m — 
2,4,6. The constant Vq = 2xoeV(ea) = 0.914ey/e, for 
Xo = 0.25. 

Table I: Madelung Matrices 





J = 2 


3 


4 




-0.4110 


0.5365 


0.347 


4^ 


-0.3466 


0.4721 


1.230 




-0.8831 


0.8831 


1.702 


K'^ — 


-0.3082 


0.6951 




K'^ — 

— 


0. 


1.082 




J^2j - 


0.1256 







The stripe phase is stable only if the_dielectric con- 
stant is large enough: recent calculations^ suggest e > 5 
is sufficient. The la*ge static dielectric constant of the 
cuprates, ^ 40 — 80Ea, is a sign of strong electron-phonon 
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coupling. This large coupling makes it difficult to accu- 
rately estimate the strength of the Coulomb interaction. 
The d.c. dielectric constant will be anisotropic and, most 
probably, dispersive, on the length scale of the stripes. 
Since interlayer contributions to screening can be impor- 
tant (e.g., apical oxygens, bilayer coupling), this is one 
parameter which could easily have a strong material de- 
pendence. 

While the above procedure should approximately cap- 
ture the long range part of the Coulomb interaction, it 
will likely overestimate the hole-hole repulsion for nearest 
neighbors. This can be thought of in terms of a correla- 
tion hole having two components. First, we are assuming 
that a hole on a given site interacts with a fractional hole 
(the average doping) on all other sites. Clearly, part of 
the hole population on the nearest neighbor sites is ac- 
tually generated by the hopping of the given hole, hence 
should not be counted in the Madelung sum. Moreover, 
there is likely to be a real correlation hole, as neighbor- 
ing charges readjust to avoid the given hole. However, 
these terms are associated with CDW formation, which 
will not be dealt with explicitly here. 



be discussed in the following subsection, which will show 
how they are modified by Coulomb interaction. 



> 




FIG. 4. Dispersion of 2,6 structure, but with structure fac- 
tors. 



IV. RESULTS 



A. Absence of Coulomb Interaction 



> 




X,Y 



S 



FIG. 3. Dispersion of 2,6 structure. 
r ^ X ^ S ^ T- dotdashed lines T ^ Y 
is along the stripes, X is across them. 



Solid lines 
S. Here, Y 



Figure || illustrates the band dispersion for a 2,6 struc- 
ture {x = 0.1875), in the absence of long-range Coulomb 
effects. The hole doping on each layer is self-consistently 
adjusted to allow for inter-row hopping processes, and the 
Fermi level is adjusted to account for the overall doping. 
In the absence of long-range Coulomb effects, the dop- 
ing is close to the nominal values. Numerical results will 
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FIG. 5. Dispersion of 2,6 structure, parallel to the stripes. 

The large number of bands is rather deceptive. It is 
equal to the number of Cu atoms in the large unit cell, 
doubled since the up and down spin bands are not degen- 
erate. There would be the same number of bands even 
if there were no stripes. But in this case, only one band 
would satisfy Bloch's theorem. This band can be deter- 
mined by looking at the structure factor - the overlap 



of the corresponding wave functions with e 



ik-r 



Simi- 



larly, when stripes are present, the same structure factor 
determines which bands will be seen by photoemission. 
This is illustrated in Fig. |, where the circles indicate a 
weight greater than 0.5, and the x 's a weight between 0.5 
and 0.1. For greater clarity, only the dispersions along 
r — > y ^ S" are shown. The dispersions along Y (par- 
allel to the stripes) and X (transverse to the stripes) 
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are shown in Figs. ^ and ^, respectively. The result- 
ing weights reveal a simple result: the envelope of the 
bands is approximately a superposition of the two lim- 
iting bands, at half filling and at optimal doping, with 
considerable fine structure associated with minigaps. 

In the presence of stripes, the dispersion should be 
quasi one-dimensional. This is clearly seen in Fig. 
where the dispersion along F — > F (F ^ X) closely re- 
sembles that along X S {Y ^ S). However, with 
the structure factors included, the dispersions are quite 
distinct. Nevertheless, the minigaps are most prominent 
in the dispersions perpendicular to the stripes, T X, 
(Fig. |), andy^5, (Fig. |). 
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mains nearly unchanged as the doping is further reduced 
(e.g., at x=0.0625). Hence, an important aspect to un- 
derstanding the strongly underdoped stripe phases will 
be to develop a good model for these limiting, two-cell 
stripes. As discussed below, there is an analogous mag- 
netic stripe beyond the percolation crossover, which can 
be modeled as a two-leg ladder. 




FIG. 7. Total dispersion from F — > S for dopings: a; = 
(a), 0.0625 (b), 0.125 (c), 0.1875 (d), and 0.25 (solid line in 
d). Data in (a) were shifted upward by 0.16eV. 



FIG. 6. Dispersion of 2,6 structure, transverse to the 
stripes. 

Fig. 1^ shows the doping dependence of the disper- 
sion, illustrating the 'projected' dispersion of the charged 
and magnetic stripes. The projected charge stripe is 
found by plotting those points where the wave function 
spends > 80% of its time on the charged layers. The 
curves show the dispersions for a series of dopings, from 
F — » X ^ 5 = (tt, tt). While the fine structure (mini- 
gaps) is strongly doping dependent, the overall dispersion 
is not, and is essentially identical to the dispersion of the 
uniform end phases. This is exactly what would be ex- 
pected for macroscopic phase separation, even though at 
crossover the charge stripe is only two cells wide. 

The figure shows that the dispersion is largely a su- 
perposition of two Fermi surfaces: one for the insulating 
magnetic stripes, one for the charged stripes: the small 
-f 's (large circles) indicate > 80% of the wave function 
is on the magnetic (charged) stripe; the small diamonds 
indicate a mixture of both. Note that there is strong 
overlap in the region of the upper Hubbard band, while 
the magnetic lower Hubbard band (LHB) remains well 
defined at all dopings, and the charged stripes fill in the 
gap as doping increases. 

It should be noted that once the charged stripes are 
reduced to two cells wide, at a;=0.125, the dispersion re- 



B. Coulomb Interaction 

Inclusion of Coulomb interaction leads to fairly modest 
changes in the dispersion. Figure 0, with no Coulomb 
effects, should be compared to Fig. ||, with moderate 
screening e = 15. Careful inspection reveals that the 
charged bands are shifted to lower energy with respect 
to the magnetic layers, so that the lower Hubbard band 
is more fully hybridized with the charged layers. With 
reduced screening (e = 5) the bottom of the charged 
band actually falls below the magnetic lower Hubbard 
band. The layers near the Fermi level remain predom- 
inantly associated with the charged layers, so we may 
still loosely speak of charged bands and magnetic bands. 
Note that in every case, the Fermi level lies within the 
minigap closest to the Van Hove singularity. This pro- 
vides a new mechanism for the opening of the pseudogap, 
as will be discussed further below. 

Even in the absence of Coulomb interaction, the car- 
rier density in a given row deviates somewhat from the 
free energy minima - here taken as x= 0, 0.25 - due 
to the finite hopping probability. For the 2,6 struc- 
ture. Figs. I^H, the magnetic layers have a;=0.025, and 
for the charged layers, moving away from the magnetic 
layer, the hole doping is 0.25, 0.24, and 0.24. Adding 
the Madelung potential raises the energy of the hole- 
doped stripes, and requires a shift of the Fermi energy 



7 



to compensate. However, since the magnetic stripes are 
gapped, this shift makes httle difference to the hole pop- 
ulation on these stripes, the corresponding layer popu- 
lations being 0.026, 0.255, 0.25, and 0.22, for e = 15, 
Fig. ^. For larger Coulomb interaction, the deviation be- 
comes greater, Fig. |[ The data display an interesting 
evolution: superimposed on a trend toward greater ho- 
mogeneity, there is also a tendency to evolve into a (2,2) 
state. This can be understood from Table I: the Coulomb 
effects are smallest for this state, since the phase separa- 
tion is restricted to the finest scale. 



.(a) 



/-(b) 




FIG. 8. Total dispersion corrected for cliarging (e — 15) 
from r ^ 5* for dopings: x = 0.0625 (a), 0.125 (b), 0.167 (c), 
0.1875 (d), and 0.25 (solid line in d). 
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0.00 



FIG. 9. Hole distribution on rows (labelled by N) of the 
2,6 structure, for e = oo (i.e., no Coulomb interaction - solid 
line), 15 (dashed line), or 5 (dotted line). 

This result is of potential relevance for LSCO: it is 
found experimentally that the incommensurability satu- 
rates near 1 /8 doping - here the crossover where the (2,2) 
phase is stable. The saturation could simply mean that 



for LSCO, the Coulomb effects are large enough that the 
system locks into the (2,2) phase for all higher dopings. 

There is a striking asymmetry about 1/8 doping: in 
the (6,2) phase. Coulomb interaction makes very little 
difference. This is because of the sharp cusp instability 
at half filling, which keeps the hole doping fixed near zero 
in the magnetic stripes, whereas the shallower potential 
minimum near optimal doping allows more substantial 
density fluctuations. 

These results will be discussed further in a later sec- 
tion. Two points are worth mentioning: first, the incom- 
mensurability saturation has so far only been observed in 
LSCO; and secondly, LSCO closely resembles the other 
cuprates in the doping range up to 1/8, but for higher 
doping, Tc saturates at a much lower value. 



C. Minigaps 

Figure ^ shows the evolution of the minigaps with dop- 
ing. A simple model provides a semiquantitative explana- 
tion of these results. The dispersion along (0, 0) (tt, 0) 
is discretized into n levels for n-Cu wide charge stripes. 
This bandwidth is A{t + 2t') ~ bSimeV. If the minibands 
are equally spaced, the average gap should be 584/ (n— 1) 
meV. Actually, the net bandwidth changes some with 
doping, so a better formula is 



bSimeV 



(9) 



= 292 (n=2), 146 (4), or 97 (6) meV, to be compared 
with average values (Fig. |8|) of 260, 147, and 94 meV, 
respectively. For the dispersion along (0, tt) — > (7r,7r), 
the same bands are present, but shifted by the dispersion 
along F, and with total bandwidth A{t — 2t'). 

For fluctuations in the stripe spacing, there will be a 
tendency to average over the various dispersions in Fig. 
^. This will tend to wash out most of the minigaps, since 
they are shifted in energy as the stripe width changes. 
However, since there is always one gap present near the 
Fermi level, this gap should survive averaging. For the 
uniform stripe phases of Fig. ^, this 'pseudogap', or dis- 
tance between the Fermi level and the nearest (tt, 0) mini- 
gap, follows the same scaling as Eq. ^, Ap = 364/ri meV. 



D. Ferromagnetic Stripes 

Figure |^ shows that, beyond the percolation crossover 
a clear remnant of the bulk VHS is visible in the striped 
phase dispersion. In Figure |l^, it can be seen that split- 
ting this VHS degeneracy produces a clear pseudogap- like 
splitting of the dispersion near (7r,0). It is this lowering 
of a large density of states that has been postulated to 
stabilize the charged stripes, and the figure clearly shows 
that the mechanism remains active even in the striped 
phase. 
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For the calculations in the figure, it was assumed that 
the antiferromagnetic phase is stable only up to a doping 
xq/S, while for larger doping a ferromagnetic instability 
wins out. The ferromagnetic dispersion is also given by 
Eqs. |-|, but with q = (0, 0) instead of Q = (tt, tt). For 
the same value of U, the equilibrium M has the form 
shown in Fig. |ll|, which was approximated by M = 0.4 — 
0.5|a;- 0.21. 




FIG. 10. Dispersion along X for a ferromagnetic instability 
on the charged stripes, for x — 0.1875 (2,6), and different 
degrees of magnetization M, as discussed in the text. 
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bilities, but one must carefully analyze secondary charac- 
teristics, since the opening of the pseudogap near (tt, 0) 
is common to a variety of instabilities. In the present 
instance, a ferromagnetic instability does not double the 
unit cell, so the ghost dispersion beyond (tt, 0) is absent, 
in contrast to experiment (see Fig. |2J below). Moreover, 
the splitting of the spin up and spin down bands should 
lead to extra structure most clearly seen (below the Fermi 
level) near F, which is not found experimentally. 



V. COMPARISON TO EXPERIMENT 

A. Photoemission in LSCO 

The doaing dependence of the photoemission spectca 
in LSCON is strikingly different from that in BSCCOEl. 
In this section, it will be shown that both spectra can 
be interpreted in terms of stripe phases, with stronger 
fluctuation effects in BSCCO. 



x = 0.1875 



FIG. 11. Doping dependence of the magnetization M for a 
ferromagnetic instability. 
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It should be noted that the doping dependence depends 
sensitively on the choice of parameters; these values are 
taken for illustrative purposes only. Figure |lo| d shows 
the dispersion of the self-consistent solution with the full 
M{x), while the other frames show a reduced M of 1/5 
(c), 1/10 (b), or (a). Since the parameters were chosen 
to have the VHS in frame (a) centered on the Fermi level, 
the pseudogap opens approximately symmetrically about 
the Fermi level. 

This should not be taken as evidence that the charged 
stripes really are ferromagnetic, only as an example of yet 
another kind of instability that is driven by the VHS. The 
figure illustrates that one can distinguish different insta- 



FIG. 12. Superposition of X and Y dispersions for, from 
bottom to top, a; = 0, 0.0625 (6,2 structure), 0.125 (2,2), and 
0.1875 (2,6). 

Since no-one has yet succeeded in making a single- 
domain sample, the photoemission should best be com- 
pared with a superposition of the X and Y dispersions. 
Figure |l2| illustrates how the pattern changes with dop- 
ing. The following features should be noted: (1) there is 
always a prominent flat band near (tt, 0), which tends to 
shift further below the Fermi level with increased under- 
doping. (2) By construction, the dispersion at half filling 
matches that found in SCOC. (3) The evolution with 
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doping is not smooth, but crosses over from the SCOC 
band near half filling to a more metallic band near the 
Fermi level at higher dopings. The overall doping is-auite 
similar to the experimental results of Ino, et al.E3 (see 
particularly their Fig. 3), confirming the suggestion that 
stripes are better developed in LSCO than in BSCCO. 

Figure |l^ corrmares the (tt, 0) photoemission peak po- 
sitions for LjSpOLj with the present calculations. Results 
for BSCCOEa are also shown; these will be discussed in 
the following subsection. In LSCO, there are two main 
features: one (x's) is near -0.6eV, with a dispersion sim- 
ilar to that in the magnetic insulator SCOC, and with a 
nearly doping-independent binding energy. The second 
feature (open circles) is a gap close to Ep with larger 
doping dependence. Qualitatively, these features are sim- 
ilar to the hump and peak features in BSCCO, but with 
larger binding energies. These two features can be cor- 
related with two prominent gap-like features in the cal- 
culations: the magnetic gap associated with the lower 
Hubbard band on the magnetic stripes, and the charge 
stripe gap, associated with the miniband closest to the 
Fermi level. The calculated gaps are larger, since the en- 
ergy scale has been chosen to agree with the magnetic gap 
in SCOC, yielding a value -1.2eV at half filling, but the 
overall doping dependences are quite similar to LSCO. 
This similarity is brought out most clearly by plotting 
the calculated gap values divided by two (diamonds and 
suns). 



B. Photoemission in BSCCO 

1. Below Tc 

The photoemission in BSCCO is strikingly different 
from the above results, yet also provides evidence for 
stripes, but of a more fluctuating form. In BSCCO, 
there is a remarkable evolution of the photoemission with 
temperature, particularly on passing through Tc. Above 
Tc, the spectra are very broad, with a single broad peak 
near tt, 0) representing the normal-state pseudogap. Be- 
low Tc, the spectra sharpen and split into two features, 
commonly referred to as a 'peak' near Ep, with a 'hump' 
at lower energies, close to the normal state pseudogap; 
between the peak and hump, there is a clear 'dip' in in- 
tensity, below the level in the normal state. Recenthi, sys- 
tematic studies-of these features in both tunnelingEj and 
photoemissiorH were presented. Most strikingly, pho- 
toemission finds these two peaks in the same direction 
of A:— space, a feature which is very suggestive of phase 
separation. 

Here, it will be assumed that the photoemission is 
dominated by stripe effects, and the main role of su- 
perconductivity is to suppress fluctuations. (The clear 
sharpening of the spectra below Tc, even in a range away 
from any gaps, is demonstrated in Ref. 70,) The analy- 
sis will be in two parts. First, the low-T spectra will be 



compared with those of LSCO. Then the role of fluctu- 
ations in producing the high-T smeared spectra will be 
discussed. 
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FIG. 13. Pseudegaps at (7r,0) in LSCOEI (x's, open cir- 
cles) and BSCCOea {+'s = hump, squares = peak) compared 
to calculated Mott gap (diamonds) and minigap (suns); the 
calculated values are reduced by a factor of two. At the high- 
est doping, the Mott gap does not show much intensity near 
(7r,0); what is plotted is energy of the corresponding disper- 
sion nearest to (7r,0). 

In p£igure ^ the low-temperature photoemission 
peaksH of BSCCO are compared to those of LSCO. The 
'peaks' (squares) are in reasonable agreement with the 
neai-Ep pseudogap in LSCO (circles), and with the cal- 
culated minigaps. On the other hand, the 'humps' (4-'s) 
are considerably closer to the Fermi level than the mag- 
netic stripe feature in LSCO (x's); compared to theory, 
the overall offset is different, but the doping dependence 
is similar. Nevertheless, identification of the hump with 
the antiferromagnetic Mott gap feature is compelling. 
LaughliiO clearly showed that the photoemission data 
evolve with doping to match the SCOC spectrum at half 
filling. A detailed fit shows that this only works if the 
dispersion at (7r/2,7r/2) touches the Fermi level - unlike 
SCOC, where there is a ^ 0.8eV gap or LSCO with a 
^ 0.3eV gap. In Appendix B it is shown that the com- 
plete doping dependence of the hump is well described 
by simply doping into the lower Hubbard band of the 
antifcrromagnet . 

There is a clear progression in the magnetic stripe sepa- 
ration from the Fermi level, from SCOC (and the present 
magnetic stripe calculation) to LSCO to BSCCO. It is 
likely that this is due to stripe fluctuation effects, since 
as the band moves closer to Ep, there will be progres- 
sively more holes on the magnetic stripes. Further evi- 
dence for this interpretation lies in the high-T BSCCO 
spectra, where the peak and hump collapse into a single 
feature, which continues to resemble the magnetic dis- 
persion (Appendix B) and is even closer to Ep than the 
hump. The conclusion that fluctuations are strongest in 
BSCCO is consistent with the fact that incommensurate 
magnetic modulations have not yet been clearly seen in 
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BSCCO. 



2. Above Tc 

In BSCCO, there is a sudden change of the photoe- 
mission spectrum at Tci a single broad feature above Tc 
spHts into a peak-dip- hump structure below T^. This is 
here assumed to be mainly a fluctuation effect: above 
Tc, both features are assumed to be present, but the line 
broadening is so large that they strongly overlap. Below 
Tc, fluctuations are greatly suppressed, and the linewidtli 
broadening F is reduced by over an order of magnitudeO, 
so a truer picture of the spectra is obtained. Supercon- 
ductivity will also renormalize the gap on the charged 
stripes (as discussed for a related model in Ref. |7|), but 
this will be a secondary effect. I— , 

In LSCO, the photoemission spectra were observedE3 
only in the superconducting state, due to surface degra- 
dation at higher temperatures. However, the stripes are 
clearly better defined in LSCO (the two gap features are 
more clearly separated) even though Tc is considerably 
lower, so it is quite possible that the split spectral peaks 
persist above Tc. 

Ref. H (Fig. 19) showed that the changes in BSCCO at 
Tc can be interpreted in terms of enhanced fluctuations in 
the normal state; a sirnilar calculation was presented by 
Chubukov and MoriO. To explain the observed BSCCO 
photoemission, it is assumed that the local hole density is 
inhomogeneous, and the photoemission can be described 
as a superposition of the spectra of different densities. 
The spectral broadening in the superconducting state is 
taken to have the normal-state value for energies greater 
than twice the superconducting gap, A^, but to have a 
(5x) smaller value at lower energies. This spectral sharp- 
ening can produce a double-peaked photoemission spec- 
trum, even if the density distribution is single peaked. A 
detailed comparison with experiment. Fig. 4a of Ref. |6^, 
will require a more detailed model of the fluctuations. 

There is an alternative class of models for the peak-dip- 
hump structure. In this picture, the dispersion is split 
because the holes strongly interact |_with some bosonic 
excitation. Abanov and Chubukovu have related the 
dip position to the resonance peak seen in neutron scat- 
tering. It is tempting to speculate that both views are 
approximately correct, and that the bosonic excitations 
originate from stripe fluctuations. 

VI. DUALITY CROSSOVER: A PERCOLATION 
QCP 

There is considerable evidgace for a quantum critical 
point (QCP) in the cupratesO. However, the QCP is at 
too high a doping to be a conventional Neel QCP, which 
should fall at a doping near x ~ 0.02, where T/v — > 0. 



Moreover, in undoped La2Cu04 and other spin-1/2 anti- 
ferromagnets, no evidence iSrfound for a high-T crossover 
to a quantum critical stateEs. Xjae QCP also lies sig- 
nificantly below optimal dopingEZl. The QCP has also 
been associated with the termination of stripe orderingL3, 
but again, the doping seems wrong: stripes persist well 
above optimal doping in Nd-substituted LSCO, as does 
the pseudogap. The evidence seems tip point to the 1/8 
anomaly playing a role: Hunt, et al.cl find a significant 
crossover, which they suggest is the end of the stripe 
regime, at x ~ 0.125, and they furtlier suggest that the 
high-field metal-insulator transitionCa takes place at this 
point, rather than at optimal doping, as reported earlier. 

The present model suggests an attractive alternative 
interpretation for the QCP, as a duality crossover of the 
striped phase. (For other patterns of nanoscale phase 
separation, such as islands, this would correspond to a 
percolation crossover, but for ID stripes, there is no true 
percolation.) At this doping, both magnetic and charged 
stripes have their minimum width, two cells. For lower 
doping, the magnetic cells widen, reducing coupling be- 
tween the charged stripes, while for higher doping, it is 
the magnetic stripes that decouple. Correspondingly, be- 
low X = 0.125, there is a charge gap, leading to high- field 
localization, while above 0.125 there is a spin gap. 

Note that this is not a conventional QCP, where there 
is an abrupt change of groundstate at T=0. For instance, 
for the 2D Hubbard model, on one side of the QCP there 
is believed to be a renormalized classical regime, with 
finite Neel order, and on the other side a quantum disor- 
dered regime, with a spin zero groundstate and a finite 
gap to the lowest triplet excitation. Instead, in the stripe 
crossover, magnetic layers persist on both sides of the du- 
ality crossover. 

While both phases persist across the duality point, 
there is a crossover in the nature of the majority phase 
(magnetic or charged), and consequently the properties 
of the minority phase are strongly modified. Because of 
this, modifications need not occur exactly at the phase 
boundary. Instead, there strong changes in properties 
may be observed at any of a series of magic dopings, cor- 
responding to the commensurate stripe phases discussed 
above. 

For instance, the localization transition in LSCO oc- 
curs not exactly at x=0.125, but closer to 0.17. This 
could be understood as evidence that two-leg charged 
stripes are always localized, but 4-leg stripes are not. 
For a uniform stripe phase, the two-leg stripes would just 
disappear at x — 2xc/3 ~ 0.167. However, the localiza- 
tion behavior does not seem to be universal: in heavily 
underdoped YBCO, the resistivity of the normal state 
satueates at low temperatures, suggestive of a metallic 
stateka. 

Again, the horizontal-vertical stripes are replaced by 
diagonal stripes at a doping 0.058, where the supercon- 
ducting transition terminates. This is close to the dop- 
ing of the (6,2) phase. (Note that the precise value of 
the magic dopings depends on xq; the pure (6,2) phase 



11 



would fall at x=0.058 if xq = 0.23.) 

Similarly, while there is not a conventional quantum 
disordered regime, a spin gap can still arise when the 
magnetic phase is the minority phase, although not nec- 
essarily at x=0.125. Indeed, in LSCO at x=0.12 ther 
is well-defined long-range (incommensurate) Neel order 
with Tm — Tc, the superconducting transition temper- 
ature. A simple interpretation of this result is that 
superconductivity is predominantly associated with the 
charged stripes, and the superconducting transition en- 
hances charge phase stiffness, reducing the fluctuations of 
the charged stripes which were suppressing magnetic or- 
der on the magnetic stripes. At higher doping, the mag- 
netic stripes are gradually spread apart, behaving more 
like two-leg ladders, and if the inter-ladder coupling be- 
comes sufficiently weak, a spin gap can open up on each 
ladder. This is discussed further in Subsection VILA. 

Note that, while a duality crossover is a generic feature 
in a stripe model based-Qn phase separation, there are al- 
ternative stripe modelS wherein the charged stripes are 
merely domain walls of antiferromagnetic domains, and 
a charge-stripe dominated regime would be meaningless. 



VII. ISOLATED MAGNETIC STRIPE 

In the stripe phase away from the duality crossover, 
the minority phase is present in the form of domain walls 
between domains of the majority phase. An important 
aspect of stripe phase theory is the development of a 
microscopic model for these domain walls. A number 
of groups have suggested a connection between magnetic 
stripes and even-leg ladders. Here, the magnetic domain 
walls in the higher-doping regime are modelled as two-leg 
ladders, which develop a spin gap as they move further 
apart, with reduced inter ladder coupling. 



straight line illustrates a modified TokuraEJ expression, 
with the doping of the planes starting ai y — 0.2, and 
varying linearly with y. The results are consistent with 
the picture that all magnetic scattering is associated with 
the magnetic stripes, and the stripe phase would termi- 
nate at an (inaccessible) doping y ~ 1.095. This would 
place the percolation crossover at ?/ ~ 0.65, close to the 
plateau regime- j-Sjince the plateau has been interpreted 
as a 1/8 effectEfUlj, this suggests that the plateau doping 
is ~ 0.125. This fixes the constant of proportionality: 
X = 0.27 {y — 0.2), so the charged stripe doping, corre- 
sponding to y — 1.095, would be ~ 0.25, in excellent 
agreement with our other estimates. At optimal doping, 
y ~ 0.925, the hole doping would be ^ 0.2_. These esti- 
mates are also consistent with Tokura, et al£3, who found 
X = 0.125 for y = 0.75, x = 0.25 for y = l,and x = .21 
for optimal doping. The inset to the fispre shows that 
Tc(x) follows the familiar parabolic formE3, 



T, 
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with Tc^max = 92X, Xm = 0.2, and Xw = 0.16. Note that 
the dip in Tc near the 60K plateau is close to a; = 1/8. 



A. Spin Gap 




In a stripe model, the magnetic neutron scattering near 
(tt, tt) should be reflective of the properties of the mag- 
netic stripes. For LSCO, the incommensurability has 
been discussed above. Fig. ^, and is further discussed in 
the following subsection. In YBCO, incommensurability 
has only been resolved at one dopingElcl, but the doping 
dependence of the peak width is consistent with a-iiimi- 
lar underlying, but unresolved incommensurabilityll3. In 
YBCO, the stripe model can also explain the doping de- 
pendence of the intensity of the magnetic neutron scat- 
tering near (tt, tt), as well as the opening of a spin gap. 

The doping dependence of the net intensity of the mag- 
netic neutron scattering should reflect the relative density 
of magneticb stripes. For YBCOg+y, the intensity was 
numerically integrated from Fig. 2 of Ref. and the 
result plotted in Fig. |l^. While the relation between y 
and hole doping x in YBCO is not completely settled, the 



FIG. 14. Intensity of magnstic scattering vs. doping for 
YBCO. Circles — data of Ref.LJ; line — expected result for 
stripe model. 

The falloff of intensity is likely to be even steeper than 
illustrated in Fig. |l^, since the data were presented only 
up to about 60meV, while at the lower doping levels, 
there is considerable intensity at higher frequencies. 

The spectrum of the excitations near (tt, tt) has a com- 
plicated evolution with doping, and below the supercon- 
ducting Tc, the intensity is suppressed Jaelow a doping 
dependent energy, called the 'spin gap'o. This gap is 
distinct from the pseudogap, and has a strikingly differ- 
ent dopiiMj .-dependence. Fig. |l^. A similar gap is seen 
in LSCOEa-Ell, but so far only near ap^p 0.15. There it is 
found that the spin gap is isotropiccZI, further evidence 
that it is distinct from superconductivity or the pseudo- 
gap. 
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The doping dependence of this spin gap in YBCO can 
be interpreted simply in terms of coupled magnetic lad- 
ders, Fig. |l5|. Below the 1/8 crossover, the magnetic 
stripe (ladder) width decreases smoothly with doping, 
while the interladder coupling is approximately constant, 
since the hole-doped stripe has fixed width. Theoret- 
ically, the spin gap is found to be (approximately) in- 
versely proportional to the ladder widthEj, so in this 
regime the spin gap scales linearly with doping. As = 
PJ/M, where J is the exchange constant, M the lad- 
der width, and /? a correction for interladder coupling, 
P ~ (1 — 4 J"/ J)|_with J" the exchange coupling between 
adjacent laddersE^. The solid line in Fig. 15 corresponds 
to J" = 0.21J. 

Above the crossover, x > xo/2 — 0.125, M is fixed at 2 
while (3 increases with doping, since J" decreases as the 
hole-doped stripes widen. Since the Cu in the hole-doped 
stripes can be magnetized, the falloff should be relatively 
slow. Details are model sensitive, but qualitatively the 
observed behavior is readily reproduced. The curve in 
Fig. [ij follows from assuming a falloff J" ~ A^~^, where 
N is the hole-doped stripe width, inset to Fig. |l^. (It 
should be noted that the falloff is sensitive to the hole- 
density xo, here takenj-as 0.25.) 

The model predicta^J that for an isolated stripe, the 
spin gap equals J/2, at least when the exchange con- 
stant is the same on all rungs and links. From Fig. |l^, 
this implies a limiting value J ~ SOmeV sX x = 0.25, 
considerably smaller than the x — Q value J — l?>QmeV . 
Such a doping dependence for J is not unexpected. For 
simplicity, however, the model assumes a constant value 
for J; this value must be taken as J = SOmeV, to suc- 
cessfully model the single stripe limit. The value is less 
critical near zero doping, where the gap is small. 
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FIG|:-45. Spin gap As vs. doping for YBCO. Circles = data 
of RefEj; line — theory, assuming solid line from inset. Inset: 
interladder exchange vs. hole-doped stripe width. 

In this model, the spin gap should already exist in 
the normal state. The striking change observed at 
can be explained as a fluctuation effect, similar to those 
seen in the BSCCO photoemission. Strong fluctuations 
at high temperatures prevent any long-range stripe or- 



der or true spin gap. The superconducting transition 
leads to three-dimensional coherence, and hence greatly 
suppresses charge and spin fluctuations in the stripes. 
Hence, a long range spin gap can open on the magnetic 
stripes below Tc- Consistent with this interpretation, it 
should be noted that the spectrum in the n prma l state in 
heavily doped YBCO has been interpr|et,edL3EZl in terms 
of a formula derived for spin-1 chainsLiJ, and hence ex- 
pected to approximately hold for spin 1/2, 2-leg ladders. 



B. Incommensurability Saturation 

Inelastic magnetic neutron diffraction finds a satura- 
tion of the incommensurabilty in LSCO at approximately 
X = 1/8. Within the present framework, there are actu- 
ally several possible explanations for the saturation. One 
was discussed above. Fig. ^: strong Coulomb interactions 
arrest the phase separation at the (2,2) stripe, and higher 
doping causes these stripes to gradually fill in. 

On this interpretation, the Coulomb effects are much 
stronger in LSCO than in YBCO, and the isolated spin 
gap regime (right-hand side of Fig. |l5|) would not exist in 
LSCO. The Coulomb effects would indeed be expected to 
be stronger in LSCO, since interlayer screening is weaker. 
In Nd-substituted LSCO, due to the LTT phase struc- 
tural distortions, stripes in alternate layers are rotated 
by 90°; a similar situation may arise in LSCO, perhaps 
due to LTT fluctuations. On the other hand, in YBCO 
the magnetic correlations have a strong c-axis modula- 
tion, suggesting that stripes in both Cu02 planes of a 
bilayer run parallel, with the charged stripes offset later- 
ally to provide stronger interlayer screening. 

However, there are other plausible explanations for in- 
commensurability saturation. Even before the discov- 
ery of stripes, it was found that LTO and LTT domains 
of fairly large size (producing distinguishiblc diffraction 
peaks) coexist near 1/8 doping in LBCO. It seems plausi- 
ble that this is associated with a stripe commensurability 
effect, similar to that found in the nickelates, and that a 
similar effect arises, at least incipiently in LSCO. In this 
case, the residual magnetic scattering would be due to 
regions that have not yet been doped beyond 1/8. In the 
more highly doped domains, the magnetic stripes would 
have a spin gap: since the ground state of a two-leg ladder 
is a spin singlet, it does not contribute to the magnetic 
scattering. A related problem has been studied by Kim, 
et al.ES, who showed that in a random mix of weakly cou- 
pled three-leg (magnetic) ladders and two-legged (spin- 
gapped, and hence non-magnetic) ladders, the magnetic 
incommensurability remains unchanged from that of the 
pure array of three-legged ladders. 

Even without commensurability effects, one would ex- 
pect 1/8 lock-in over a finite doping range, when the 
2,2 stripes coexist with 2,4 stripes, which have a well 
defined spin gap. In this case, the magnetic incommen- 
surability should be fixed at that for 1/8 throughout the 
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coexistence regime, but should disappear when a com- 
mensurate 2,4 phase is stable, at a; = 2xn73 ~ 0.17. In 
LSCO, the 1/8 stripes are actually foundElfi to persist up 
to X ~ 0.25. Hence the need to postulate lock-in cflccts 
at 1/8 doping, the exact analog of the stability of the 
X = 1/3 and 1/2 striped phases in nickelates. In this 
case, the heavily doped phases would have no magnetic 
scattering, while the 1/8 stripes would have a scatter- 
ing of fixed incommensurability, but decreasing intensity 
and increasing width, as the stripe domains shrink in 
size. The special stability of the 1 /8 phase may be asso- 
ciated with the finite residual exchange coupling across 
the two-cell-wide charged stripes, which is responsible for 
the antiphase boundaries, and which may be lost in wider 
charged stripes, or with the reduced Coulomb energy. 

At this stage, there is not enough information to judge 
between the two models for incommensurability satura- 
tion. The former, strong Coulomb effect, has the ad- 
vantage that it could simultaneously explain why Tc in 
LSCO is so low - the local hole density is forced away 
from optimal. However, there is considerable evidence 
that stripe phase order is better developed in LSCO than 
in other cuprates, and this could provide reason enough 
for a lower Tc. 

In many ways, optimally doped LSCO resembles an 
underdoped YBCO. We have here suggested that this is 
because stripes and pseudogaps in both materials persist 
up to a; ~ 0.25, whereas Tc-iis optimal near 0.16 in LSCO, 
0.2 in YBCO. Sato, et alz3 have recently provided addi- 
tional evidence that the pseudogap opens well above Tc 
in optimally doped LSCO. 

VIII. DISCUSSION 

A. Improvements for the Stripe Model 

The present model of the stripe phase provides a sig- 
nificant advance over earlier calculations. In these early 
calculations, an external stripe potential was imposed, 
which could be either periodic or random, and the rear- 
rangement of holes was studied. It was found that the 
stripes produced minigaps, but that an average, smeared 
dispersion and Fermi surface could still be defined. How- 
ever, there was no sign of the split dospersion (magnetic 
vs charged stripe) found in photoemission experiments. 
The present calculation replaces the external potential 
with a given, doping dependent potential, and studies 
how holes redistribute in the presence of a competition 
between that potential, which favors phase separation, 
and Coulomb repulsion, which favors a uniform density 
distribution. The striking result is that the resulting dis- 
persion resembles a weighted superposition of the disper- 
sions of the two end phases, with the addition of some su- 
per lattice minigaps. This is very encouraging, in provid- 
ing an explanation for the photoemission results. More- 
over, it shows that the mechanisms responsible for the 



special stability of the end phases can continue to oper- 
ate on these nanoscopic length scales. 

A number of improvements still need to be made in 
the model. The next step would be to make the calcula- 
tion fully self-consistent, by eliminating any assumed po- 
tential, and directly calculating and minimizing the free 
energy of the striped phase. Since the dispersion is not 
greatly changed, it is unlikely that this additional step 
will greatly modify the present results. The most likely 
change would be that the densities could adjust slightly 
to take advantage of the minigaps, better centering them 
at the Fermi level. This could lead to a more system- 
atic growth of the pseudogap with underdoping, since the 
minigaps are associated with the charged stripes, and get 
larger as these stripes get narrower. The idea that the 
pseudogap is associated with stripe minigaps has been 
proposed previouslyl£j; the present calculation provides a 
systematic doping dependence and a connection with the 
VHS. 

A complete understanding of the stripe phase, partic- 
ularly in BSCCO, lies in the correct inclusion of fluctua- 
tion effects. These effects can broadly be separated into 
two categories, depending on whether the fluctuation pre- 
serves the local density distribution or not. In the former 
category fall fluctuations in the local stripe spacing, ei- 
ther static or slowly varying in time, and long-wavelength 
bending of the stripes. It is likely that the energy dis- 
persion is a fairly localized function in space, and that 
these fluctuations can be calculated as weighted averages 
over the present solutions. In this case, the dispersion 
would still be a superposition of the two end phases, and 
the main effect of the fluctuations would be to smear out 
the minigaps. Since there is always a gap near the Fermi 
level, a residual pseudogap should survive. Moreover, 
since the split-off LHB is well defined, particularly at 
lower doping, it should persist at a distinguishible feature 
after averaging. This would resemble the photoemission 
in LSCO. 

The second class of fluctuations involves fluctuations 
which are fast enough, or disordered on a sufficiently 
short-wavelength scale so that the local density does not 
lie near the two potential minima. These fluctuations 
act to wipe out the stripe fluctuations on a local level, 
and the question is, can they describe the experimen- 
tal results in BSCCO, where, for T > Tc, the two valence 
bands appear to collapse into a single reconstructed band. 
This is a plausible result: as a line of holes fluctuates back 
and forth in an antiferromagnetic background, the back- 
ground will have to adjust to some time-averaged hole 
density. The theoretical problem is how to properly in- 
clude this averaging: it is a question of how the system 
responds locally on different time scales. 
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B. Four Types of Gap 

It should be noted that the present model has four 
distinct types of gap in the electronic spectra. First, there 
is the Mott-Hubbard gap in the magnetic stripes, e.g.. 
Fig. I^a. Then there is the Van Hove gap on the charged 
stripes. Fig. |lo|. It is postulated that the opening of these 
two gaps provides the stabilization energy of the stripe 
phase, as in Fig. |^. However, in addition there is a third 
type of gap, arising from the new periodicity associated 
with the stripe phase. This leads to minigaps in the 
dispersion across the stripes, clearly visible in Fig. ^-d. 
Finally, there is the superconducting gap. 

Many earlier discussions of the pseudogap relied on 
the photoemission studies of BSCCO in tha—uarmal 
state. |Xhe recent systematic studies in BSCC0E3O and 
LSCOEa suggest that the low-temperature data are more 
representative of the stripe phase. In this case. Fig. |l^ 
shows that there are two distinct features, and the Mott- 
Hubbard gap is associated with the 'hump' feature in 
BSCCO; the 'peak' feature is then plausibly associated 
with the charged stripes. The question then is which of 
the remaining three theoretical gap features accounts for 
the 'peak'. 

The doping dependence of the peak feature is not con- 
sistent with a conventional superconducting gap: the 
transition temperature decreases while the peak en- 
ergy Ap increases with increasing underdoping. It has 
been suggestedEj that Ap is superconductivity-related, 
but associated with preformed pairs. However, in the 
case of a multicomponent order parameter, it has been 
shown that the total gap at (tt, 0) is a vector sum of the 
individual, components. In BSCCO, a careful analysis of 
the dataEj (involving the detailed doping dependence of 
the 'dip' feature) suggests that the gap is best interpreted 
as a superposition of a superconducting component and 
a non-superconducting component, with the latter dom- 
inating in the strongly underdoped limit. Accepting this 
analysis, the question becomes, what is the origin of the 
non-superconducting component on the charged stripes? 

The experimental observation that this peak increases 
as the hole doping decreases cannot be easily explained in 
terms of conventional Van Hove gaps, which are largest 
at xq, Figs. 1^ anii] O- A model involving Van Hove pin- 
ning, introduced^ to explain the pseudogap crossover in 
BSCCO, relied on the normal-state data, which appeared 
to show a smooth evolution into the magnetic stripe at 
half filling. The model cannot explain the evolution of a 
distinct charge-stripe gap. 

On the other hand, the minigaps seen in Fig. ^ do 
on average grow larger as the doping is reduced, satu- 
rating near the 1/8 crossover. In the present, non-self- 
consistent calculation, these minigaps are not centered 
on the Fermi level. However, it is clear that having Ep 
fall in a minigap would enhance the energy lowering, sta- 
bilizing these gapped phases and improving agreement 
with experiment. It should be noted that the greatest 



energy lowering would correspond to having Ep centered 
on the minigap closest to the VHS, since that is where 
the dos is highest. Fig. 16, This would then be a form of 



stripe-induced Van Hove splitting. 
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FIG. 16. Density of states for dopings: a; = (a), 0.0625 
(b), 0.0833 (c), 0.125 (d), 0.167 (e), 0.1875 (f), and 0.25 (g) 
assuming dielectric constant e — 15. (The a:: = data are 
shifted up by O.lGeV.) Solid (dashed) line — partial density 
of states for magnetic (charged) stripes. 
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FIG. 17. Blowup of density of states near Ep, for same 
dopings as in Figure 0. (There is no frame (a), since the dos 
vanishes in this energy regime.) 

Note that the largest minigap is associated with this 
VHS splitting, even in the lowest doping situation, (6,2). 
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Figure shows a blowup of the dos near the Fermi 
level. It can be seen that there is a gap (dos = 0) for 
X < 0.125), centered at an energy ~ O.lSeT^ above Ef. 
For larger x a pseudogap (dos > 0) persists, shifting to- 
ward Ep with increasing doping. Since the problem is 
two dimensional and the stripes are fluctuating, this gap 
will be spread out into a pseudogap. Nevertheless, the 
low dos means that localization effects should be present. 
Since the (pseudo)gap remains large at the 1/8 crossover, 
there should be a delocalization transition at a somewhat 
higher doping. In the present calculation, this feature 
falls above the Fermi level; however, this could change in 
a fully self consistent calculation. 

Finally, one should note the duality between the charge 
stripe minigap, which grows as hole doping is reduced, 
Fig. 13, and the magnetic stripe spin gap, which increases 
with increasing hole doping. Fig. |l^. 



C. Fermi Surface and Remnant Fermi Surface 



fall to half at an energy away from the Fermi level (N.B., 
the rFs is not a surface of constant energy). Indeed, a 
rFs was found for the insulating CCOC. We have shown 
that in this case the rFs (the locus of points where the 
coherence factor equals one half) maps out the superlat- 
tice zone boundary. (Since the model does not include 
fluctuations, there is only one band below Ep, and hence 
no photoemission distribution to integrate over.) 

In the case of a stripe array, there are several complica- 
tions. First, there are several subbands, and one will get 
different results depending on whether one calculates an 
rFs for each subband, or a single rFs for the whole valence 
band. The structure factor provides an additional com- 
plication, since the intensity is almost never the full pos- 
sible value. Nevertheless, for simplicity. Fig. |l^ plots the 
locus of points where the net spectral function equals 1/2 
- actually, falls within a range 0.48-0.52. Comparison of 
Figures |l^ and |l^ shows that the true Fermi surface and 
the rFs are quite distinct features, and that the rFs tends 
to follow the superlattice Brillouin zone boundaries. 



Figure shows the Fermi surfaces corresponding to 
the same dopings as in Fig. ^ As expected from the 'pro- 
jected' dispersions in that Figure, in all cases the wave 
functions are > 90% associated with the charged stripes. 
Note the important role of the structure factors: while 
there is a well-defined superlattice in each case, and hence 
each Fermi surface segment is periodically repeated, the 
weight is highly nonuniformly distributed, being concen- 
trated predominantly near the limiting x=0.25 Fermi sur- 
face (solid IhiCpin the figures). 

Feng, et al.EZI recently presented photoemission evi- 
dence for the presence of two Fermi surfaces in BSCCO, 
seen at different incident photon energies. This result is 
disputedEa, but it is intriguing that the 'new' Fermi sur- 
face has large flat sections, very similar to the 1/8 stripe 
pattern, Fig. ^8|d. 
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FIG. 18. Fermi Surfaces for dopings: x = 0.0625 (a), 0.125 
(b), and 0.1875 (c), and e = 15. Solid line in each is the Fermi 
surface for x=0.25. 

Ronning, et al.S introduced an alternative, well- 
defined energy surface, which they refer to as a 'remnant 
Fermi surface' (rFs). This is the locus of points where the 
integrated photoemission intensity, taken as proportional 
to n{k), falls to one half its maximal value. While the 
intensity does fall to half at the Fermi level, it can also 
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FIG. 19. Remnant Fermi Surfaces for same dopings as in 
Fig. |l|: X = 0.0625 (a), 0.125 (b), and 0.1875 (c). 

D. Staging 

A test case for phase separation models of stripe 
formation comes from measurements on oxygen-doped 
La2Cu04-f5. In these materials, the interstitial oxygens 
are highly mobile, so the tendency of holes to phase 
separate leads to macroscopic phase separation. As in 
graphite intercalation compounds (GIC), this is accom- 
plished by stagingO: the interstitial oxygens form dense 
layers, with subsequent oxygen layers separated by n 
layers of CU.O2 (n is the stage number). The puzzling 
observationEl is that despite the macroscopic phase sepa- 
ration, stripes are still found in these compounds. I sug- 
gest the following resolution, also based on an analogy 
with GIC. The stripes were observed in a predominantly 
stage 4 compound, with 4 Cu02 layers per oxygen inter- 
stitial layer. In GIC, a stage 4 would have a highly inho- 
mogeneous charge distribution, with most of the charge 
on the bounding layers (in this case, the Cu02 layers ad- 
jacent to the interstitial O's), and significantly less on the 
interior layers, away from the O's. Hence, in the cuprate I 
propose a similar effect: fully charged x ~ 0.25 bounding 
layers, and lighter charged '--^ 1/8 interior layers. 
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Thus, the staging is driven by the tendency to phase 
separate, producing ideal hole-doped bounding layers, 
but there are still stripes due to doping of the inte- 
rior layer. This model resolves two additional problems. 
First, the sample was a mixture of stages 2 and 4. If the 
doped holes were uniformly distributed on all Cu02 lay- 
ers, then there would be two independent stripe patterns, 
since the stage 2 Cu02 planes would have nearly twice as 
many holes per layer. Within the present model, most of 
the charge goes into the bounding layer, so the layers of 
stage 2 would have comparable doping to the bounding 
layers of stage 4; moreover, since both types of layers are 
nearly fully doped, neither would show magnetic scatter- 
ing. 

Finally, even though static stripe order was observed, 
there was no reduction in Tc - indeed, Tc=42K is higher 
than can be obtained with Sr doping. This can most 
easily be understood if the superconductivity and the 
stripe order are in different layers. 



E. Comparisons with Slave Boson Results 

It is instructive to compare the present results with 
earlier slave boson calculations. In the simplest version, 
there is no magnetic coupling, J = 0, and the band struc- 
ture near the Mott transition is highly anomalous. There 
is a single band, but as the doping approaches half fill- 
ing, X — > 0"*", the bandwidth vanishes, with both t and t' 
renormalized to zero. In the three-band model, even af- 
ter setting J7 — !■ oo, there is still a charge transfer energy, 
A In this case, it is also possible to approach half filling 
from below, a; — > 0~; the same bandwidth collapse oc- 
curs, but at a different energy, E~ , with — being 
the (renormalized) charge transfer energy. 

Is there any way to reconcile the present results with 
slave boson theory? I suggest the following possibility. 
When a hole is doped into the Mott insulator, there is 
phase separation, and locally the dispersion is restored: 
t tQ. At a doping x, a fraction x/xq of the electrons 
have hopping ~ t, the rest ~ 0. But in the mean-field 
slave boson calculation the effect of the hole is uniformly 
spread out over the entire lattice, leading to an effective 
t xtQ. This is just what is found in the present stripe 
calculation. As the material is doped, the magnetic band 
persists with little change, while a new band appears, 
characteristic of the hole-doped stripes, with full band- 
width, t ^ to (neglecting super lattice gaps), but with 
relative intensity proportional to x, Fig. ^ If this in- 
terpretation is correct, it suggests that the slave boson 
calculation may underestimate the tendency for phase 
separation. 



F. Comparisons with Other Calculations 



A number of calculationalHl 



interpret the mag- 
netic stripe dispersion in terms of a flux phase rather 
than an antiferromagnet. This is unsatisfactory for two 
reasons: first, there is clear evidence that the mag- 
netic stripes have predominantly antiferromagnetic cor- 
relations; and secondly, the flux phase has zero gap at 
(7r/2,7r/2), making it difficult to explain the dispersion 
in SCOC, which has a gap > 0.8eV. The present cal- 
culation is based on an antiferromagnetic model, which 
corrects both of these defects. 

The dispersion inj-the presence of stripes has-also been 
calculated by SEKEJ and by Seibold, et aLEll. These 
calculations were intended simply to show the effect of 
charge (or spin) modulations on the quasiparticle disper- 
sion and involve highly simplified models for stripes. SEK 
introduced a background stripe potential, Va-{R), which 
fixed the orientation and periodicity of the stripe array, 
and induced a modulation in the quasiparticle charge and 
spin densities. Seibold, et al. concentrate on the charge 
modulation, and introduce a CDW model, again with 
externally imposed periodicity. 

While these calculations demonstrated how the mod- 
ulation introduces pseudogaps into the dispersion, they 
are not expected to give a very good picture of stripes 
arising from an underlying phase separation instability. 
A CDW model involves a fixed periodicity and a variable 
density, whereas phase separation should produce an ap- 
proximately fixed density and smoothly varying period- 
icity. Thus, in the present calculations, it is straightfor- 
ward to model the full doping dependence of the stripe 
phase; in the earlier calculations, the doping dependence 
of the periodicity and potential must be supplied em- 
pirically, so the models can make no predictions. Thus, 
neither calculation showed any sign of the separate mag- 
netic and charge bands. 

One should note a distinction between the way the idea 
of CDW's is applied in these calculations and tlift-usage 
of the present paper. In the earlier calculationscJ'EJ, the 
charge modulation of the stripes is the CDW. In contrast, 
the present paper envisages CDW's (or related phonon 
anomalies, such as the LTT phase) as existing on the 
charged stripes, leading to two inequivalent sublattices - 
much as antiferromagnetism leads to two sublattices on 
the magnetic stripes. 



In a recent paper, Pryadko, et al.trfl presented a de- 
tailed criticism of the early Hartree-Fock calculations of 
stripes. In these calculations, the stripes are claimed to 
arise from Fermi surface nesting, so some of the criti- 
cisms apply to a larger class of models. Here, I would 
like to take issue with some of the statements. First, 
Pryadko, et al. state "the density of holes along a stripe 
varies, continuously as a function of x" , citing Yamada, 
et alfl. While this is mainly associated with the satura- 
tion of the incommensurability above x—l/8, they note 
that the data at lower doping also show some curvature, 
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suggestive of a 'slight variation' of the doping. This is 
mainly associated with the cutoff of incommensurability 
below X ~ 0.06, but that data have now been superseded 
by better samples: the stripes are present, but rotated 
by 45°, and the incommensurability falls close to the Ya- 
mada (straight) line. As discussed above, both a shght 
doping-dependence below x = 1/8, and a stronger satu- 
ration above that doping can be understood in a model 
with a prefered charged stripe density, in the presence of 
either long-range Coulomb repulsion (which can produce 
a real saturation of the incommensurability) or commen- 
surability pinning (which produces an apparent satura- 
tion). 

Pryadko, et al. further state, "In the LSCO family, 
there is simply no vestige of a quasiparticle in the re- 
gion of momentum space where the nested Fermi surface 
is supposed to occur." This statement raises several is- 
sues. First, one must carefully distinguish conventional 
Fermi surface nesting from,-\ian Hove nesting. For the 
conventional nesting modelE2^ they cite, the nesting is 
associated with flat sections ,€>£, Fermi surface; in LSCO 



these flat sections arise nearEHil (7r/2, 7r/2), where there 



are well-defined quasiparticles. Conversely, the fact that 
these quasiparticles persist in the striped phase shows 
that they are not involved in the nesting process. 

In contrast, in Van Hove nesting the important quasi- 
particles are those near (tt, 0). In all the cuprates, these 
quasiparticles are found to be shifted below the Fermi 
level, and indeed a comparison of photoemission and tun- 
neling studies reveals that the position of these quasiparti- 
cles defines the pseudogap — exactly as required by a Van 
Hove nesting theory. Moreover, even in the SEK calcu- 
lation it was clear that these Van Hove quasiparticles, 
which constitute the 'flat bands' remain well defined in 
the calculated striped phase. This is further clear from, 
e.g.. Fig. 1^ above. True the quasiparticles are greatly 
smeared out, but this is a direct consequence of the for- 
mation of stripes, and cannot be taken as an argument 
against this mechanism of strip&,formation. 



Finally, Stojkovic and PinesliH3 claimed to rule out the 



VHS in the physics of the cuprates. They employed pa- 
rameters for which the VHS falls at very high doping, 
X = 0.55, and concluded "the presence of the van Hove 
singularities near the Fermi surface plays only aujparginal 
role" . However, their more recent calculationsa employ 
a revised parameter set (i.e., t' = —0.25^), which places 
the VHS at a; = 0.22, close to the present xq — 0.25. 



IX. CONCLUSIONS 

This has been a long manuscript, which presents a 
coherent view of the stripe phases in the cuprates. A 
number of principal results of the calculations are here 
summarized. Most of the results are generic^ and would 
be expected in any model where the stripes result from 
two-phase coexistence, while a few are specific to a Van 



Hove scenario. 

(1) This is the first calculation of photoemission in the 
stripe phase involving strong stripe correlations, i.e., pre- 
ferred hole dopings with independent characteristic dis- 
persions, as opposed to a single dispersion with sinusoidal 
(CDW/SDW) modulations. 

(2) This allows a study of the evolution of the disper- 
sion as a function of hole doping. 

(3) It is found that, even at this nanoscale level, the 
dispersion can be characterized as a superposition of two 
components, leading to a picture of magnetic stripe bands 
and hole-doped ('charged') stripe bands. This allows a 
natural interpretation of the photoemission spectra in 
LSCO, and suggests a unified picture with BSCCO and 
SCOC. 

(4) The calculations suggest that an important role 
of the superconducting transition is to freeze out fluc- 
tuations of the stripes. This freezeout manifests itself 
in three ways: (a) the electron-electron scattering rate 
drops by several orders of magnitude below TcO; (b) the 
photoemission dispersion splits in BSCCO into a charac- 
teristic peak-dip-hump structure; (c) the (tt, tt) magnetic 
neutron scattering in YBCO sharpens below Tc, reveal- 
ing a characteristic spin gap. 

(5) The doping dependence naturally leads to a picture 
of a series of quantum critical points (QCP's) or magic 
dopings, at which the stripe pattern is commensurate 
with the crystalline lattice. The most prominent one is 
the famous 1/8 effect, but the metal-insulator transition 
in LSCO and the onset of superconductivity are close to 
two other magic numbers. 

(6) The percolation crossover at 1/8 doping provides 
a simple model of the spin gap in YBCO, showing that 
a two-leg ladder provides a good model for an isolated 
magnetic stripe. 

(7) As a result of point (3), the model has a natural 
VHS pinning to the Fermi level: if the VHS is at the 
Fermi level in the charged stripe end phase (as it must be, 
if this phase is stabilized by Van Hove nesting) , then the 
VHS remains close to Ep over the entire doping range. 

(8) This provides a new explanation of the pseudogap: 
stripe-induced Van Hove splitting. 

(9) More speculatively, since superconductivity in 
YBCO is strongest well beyond the percolation crossover 
(1/8 effect), superconductivity seems to be a property pre- 
dominantly of the charged stripes. 

There are a number of advantages of the present model 
of fractionally-doped stripes. First, if the stripes are sta- 
bilized by CDW formation, then there is an important 
continuity between stripes in the cuprates, and those in 
the nickelates and manganites. Such continuity is lost 
in the SO (5) model, where the charged stripes are stabi- 
lized by superconductivity. Moreover, a connection with 
CDW's would naturally explain the experimental obser- 
vation that the stripe phases are dominated by charge 
order rather than spin order, a result difhcult to under- 
stand in a pure Hubbard or tJ model. 
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more clearly seen in plots of the rung-averaged hole den- 
sity, Fig. gy. The data (open circles) fall very close to the 
form expected for a phase separation model (solid lines), 
with the same average densities as at 1 /8 filling, but now 
the magnetic stripes retain their minimum width, while 
the hole-doped stripes get wider. 



APPENDIX A: WHITE-SCALAPINO STRIPES 

It remains an open question whether phase separa- 
tion is a generic feature of the Hubbard and tJ mod- 
els, or arises only in a restricted parameter domain. 
In the tJ model, the recent density-matrix renormaliza- 
tion (DMRG) calculations of White and Scalapino 

(WS)cJ'EJ find clear evidence for stripes, but even these 
results are controversial: some groups, Sad that stripes 
are metastable within th e tJ mode£2j^E£] (see the dis- 
cussion in Refs. 108,109), others that—the model should 
display macroscopic phase separationtij. Moreover, in- 
clusion of realistic values, of./' into the model further re- 
duces the stripe stabilityllHllli3. Nevertheless, they rep- 
resent at worst a low-lying excited state, and hence may 
become the true ground state in the presence of addi- 
tional interactions (e.g., electron-phonon coupling). 

It is interesting to note that many of the DMRG results 
can be simply understood in terms of a phase separation 
model. This hypothesis can explain fcuir quantitative re- 
sults of WS: (1) At 1/8 doping, WSEI find hole-doped 
stripes which are 2 Cu atoms wide, with doping ^ 0.25 
holes per Cu, which they interpret as a bond-centered do- 
main wall with 0.5 holes per unit length. Increasing the 
doping, (2) they find a remarkable transition to stripes 
with 1 hole per unit lengthEll, (3) which coexist with the 
0.5 hole stripes for 0.125< x < 0.17, and finally, (4) the 
stripe phase disappears at x=0.3. 

All of the above features can be simply understood 
by assuming that the hole-doped stripes have a density 
X* ~ 0.25 hole/Cu, and that there is a preference for 
both magnetic and charged stripes to have even width. 
Hence, the narrowest possible magnetic stripes are two 
cells wide. 

Given these assumptions, the WS data can be ex- 
plained as follows. In the low doping regime, hole doped 
stripes form with the minimal width of two cells, and 0.25 
holes per Cu, for a net of 0.5 holes per unit length. This 
continues up to 1/8 doping, at which point both magnetic 
and hole-doped stripes are two cells wide. Doping holes 
beyond 1/8 doping can be accomplished by increasing 
the width of the hole-doped stripes. A stripe 4 Cu wide 
has 4x0.25 = 1 hole per unit length. The sample will be 
a mix of 2-Cu and 4-Cu stripes until x — 2x* /3 ^ 0.167, 
and when x = x* ~0.25, the whole sample is filled with 
the hole-doped phase, and the stripe phase terminates. 
This last number is lower than the simulation, x = 0.3, 
and may be evidence that x* is weakly doping dependent. 

The WS simulations actually display this iMrease in 
stripe width above 1/8 doping. Fig. 4a of Refo. This is 



0.25 
0.20 h 
0.15 

- 1 

0.10 " 

0.05 

0.00 




_J i_ 

/ \ 
i t 
I \ 



10 15 20 
row 



25 



FIG. 20. Rung-averaged hole density (open circles)0, in- 
terpreted as constant hole density domains of variable width. 
Dashed lines — guides to the eye; solid lines = phase separa- 
tion model. 

It should be stressed that the mechanism responsible 
for stabilizing the WS stripes is unknown, but clearly 
does not involve the VHS, being weakened by a finite t'. 
Even should the WS stripes prove to be the ground state, 
it may still be that they are too weakly bound to explain 
experimental observations. The debate on their stability 
suggests that non-striped phases lie very^lose in energy 
for instance, Hellberg and ManousakisE£3 find energy 



differences of the order of ImeV. But it is found that the 
stripe phase separation collapses when the temperature 
becomjes-.comparablc to the energy barrier between the 
phaseiia. Fig. |2l|. 
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FIG. 21. Striped phase melting transition at x=0.04, 
e = 23: hole distribution function at temperatures (reading 
from bottom to top, at fixed hole occupancy 30a; = 2) fcsT 
= 10, 0.1, 30, 60, 100, 400, 200 meV. Inset: Assumed free 
energy vs. doping. 
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It should be noted that Figure 2 of Ref. 112 was printed 
incorrectly in the original publication. The correct figure 
is included here, as Fig. |2l[ The calculations are de- 
scribed in the original publication. Note that, at low 
temperatures, the particle distribution is two peaked, 
corresponding to the two kinds of stripes, while at high 
temperatures the distribution collapses into a monotonic 
limit. The crossover temperature scales with the free en- 
ergy barrier between the two phases. 



APPENDIX B: MORE ON THE SDW 
DISPERSION 

As will be shown below, the SDW model gives a sur- 
prisingly good description of the photoemission disper- 
sion, not only in the insulating limit, but also in doped 
BSCCO. However, if M is interpreted as the long-range 
Neel order parameter, there are problems with the dop- 
ing dependence. Fig. |l|, and the temperature dependence. 
Fig. On the other hand, since this is a mean-field cal- 
culation, it is more plausible to interpret M as a measure 
of short-range antiferromagnetic order, and the transition 
as describing the splitting of the upper and lower Hub- 
bard bands. To stress this point, the transition is labelled 
as Tfj. Figure ^ shows the temperature dependence of 
the magnetization, while the inset shows the dependence 
of Th on U /t. These results are consistent jKith the three- 
dimensional calculations of Dichtel, et al.EJ. 



tent with the pseudogap : 



0.4 - 



0.3 



M 



0.2 



0.1 



0.0 



2000 







'1 


'"I""!""!""; 


































/ 






,1 





2 4 6 8 10 

u/t 

I .... I .... I .... I 



1000 2000 3000 4000 5000 

T 

FIG. 22. Temperature dependence of magnetization in 
mean-field model. Inset: 17-dependence of Th- Dashed line: 
T = U/i. 

Interpreting the doping dependence of M as a mea- 
sure of the renormalization of the splitting into upper and 
lower Hubbard bands, the mean- field calculations are ia 
good agreement with exact diagOBalization calculationsES 
and experiments on the cupratesES. 

Figures 23-26] show how the dispersion changes with 
doping. For all curves, the same band parameters, given 
below Eq. U, are used, with M{x) found self consis- 
tently. Surprisingly, the results are qualitatively consis- 



e.underdoped cuprates (x 's, 
diamonds, and squares even though no striped 
phases are assumed. However, the predicted hole dop- 
ing is anomalously large. 
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FIG. 23. Dispersion of the antiferromagnetic insulator in 
mean-field model. Open circles = data of Ref.EZI. Solid lines: 
coherence factor > 0.8; dashed lines: 0.8 > coherence factor 
> 0.2: dotted lines: coherence factor < 0.2. 
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FIG. 24. Dispersion of the doped antiferropiftgnet 



X = 0.169, in mean-field model, x's = data of Ref.lii3 (un- 
derdoped, Tc = 67K); diamonds (underdoped, Tc |=j 52K) 
and squares (overdoped, Tc — 72K) = data of Ref.l-3; solid 
lines: coherence factor > 0.8; dashed lines: 0.8 > coherence 
factor > 0.2: dotted lines: coherence factor < 0.2. Horizontal 
dashed line = Fermi level. 

In the undoped material, the model has a direct (in- 
direct) Mott-Hubbard gap of 1.63eV (1.27eV), antLm- 
produces the dispersion found in the oxycloridesEZlEj, 
Fig. l23l Here X = (tt, 0), 5* = (tt, tt), and S = S/2. The 



spectral weight is proportional to the coherence factor 
C± = i(l±^), (Bl) 



with the subscript -f- (— ) referring to the upper (lower) 
Hubbard band. From Fig. it can be seen that the 
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coherence factor is < 0.5 all along the X ^ S branch, 
where no quasiparticle peaks were reported. The weight 
along the line X ^ S is uniformly l/2|-characteristic of 
the antiferromagnctic zone boundaryc3lli3. 



shall, et aLliiJ are plotti; 
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FIG. 25. Dispersion of the doped antiferromagnet, 
X = 0.356, in mean-field model. Diamonds = data of Ref.tJ 
(overdoped, Tc = 85-?^); solid lines: coherence factor > 0.8; 
dashed lines: 0.8 > coherence factor > 0.2: dotted lines: co- 
herence factor < 0.2. Horizontal dashed line = Fermi level. 



; the normal state pseudogap, 
plot the hump feature in the 



while Campuzano, et al.l 
superconducting state.] 

For these parameters, the transition to the paramag- 
netic phase is first order; for other parameter choices, a 
second-order transition is found. The first order transi- 
tion is a topological transition, arising when the Fermi 
level crosses the band dispersion near X (Fig. |2^). This 
is rather striking, since topological transitions are typ- 
ically rather weak - of order 2.5. A similar result 
found in the Hubbard model {t' — 0) by Guinea, et al.l 

Thus, the dispersion of the normal state pseudogap, as 
seen in photoemission, can be well described by a uni- 
form doping of the magnetic phase, in the absence of 
stripes. This has been npted by a number o£,groups, in- 
cluding Schmalian, et al.til, and Misra, et al.E3. However, 
this cannot be the complete interpretation. Tunneling 
studies clearly demonstrate that the pseudogap consists 
of two components, split approximately symmetrically 
about the Fermi level. The uniform phase dispersion de- 
scribes the feature below the Fermi level, but the model 
contains no corresponding feature above the Fermi level. 
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FIG. 26. Dispersion of the doped antiferromagnet, 
X = 0.392, in mean-field model. Solid lines: coherence factor 
> 0.8; dashed lines: 0.8 > coherence factor > 0.2: dotted 
lines: coherence factor < 0.2. Horizontal dashed line = Fermi 
level. 

In studies of the Fermi surfaces of the hole-doped 
cuprates, there is an important question as to whether 
the Fermi surfaces are large, centered on F, or small 
hole pockets, centered at S = (7r/2,7r/2), and equivalent 
points. In dispersion studies, the hole pockets are associ- 
ated with retrograde dispersion branches of low intensity 
(ghosts), particularly along X ^ S and S ^ S. Note 
that this negative curvature along the branch X -jr S. is 
clearly seen for undcrdoped samples with < 75X113^3, 
Fig. 24, [Parenthetically, the striking difference in doping 
for the two data sets in this figure arises because Mar- 
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